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The  VIPER  Conductivity  Model 

INTRODUCTION 

In  this  paper  we  briefly  discuss  the  conductivity  models  used  in  the 
electron  beam  propagation  code  VIPER^ .  This  is  a  linearized  code  built  to 
study  the  resistive  hose  Instability  in  an  electron  beam  propagating  in  air. 
The  beam  Ionizes  the  air  creating  a  plasma  whose  conductivity  modifies  the 
beam  generated  fields.  VIPER  uses  the  ultrarelativist ic  frozen  approximation, 
replacing  z  and  t  by  z  (the  position  in  the  lab  frame)  and  z  •  ct-z  (the 
distance  from  the  beam  head).  In  a  typical  run  there  may  be  200  beam  slices, 
each  slice  having  50  radial  points  at  which  the  conductivity  must  be 
calculated.  These  slices  move  «  400  steps  in  space.  In  addition,  the 
perturbed  or  dipole  conductivity  must  be  found  at  all  the  points. 

A  complete  calculation  of  the  conductivity  is  extremely  difficult  . 

Dozens  of  species  of  neutrals,  ions  and  their  excited  states  are  created. 

They  react  with  the  beam,  the  plasma  electrons,  and  themselves.  In  many  cases 
the  reaction  rates  are  not  known  or  there  is  inconsistent  experimental 
t  evidence  for  them.  The  plasma  electrons  are  mobile  and  are  affected  by  the 

electric  and  magnetic  fields  of  the  beam.  In  addition,  chemical  reaction 
systems  are  notoriously  stiff  numerically.  Computational  models  which  offer  a 
comprehensive  treatment  of  air  chemistry  and  conductivity  generation,  such  as 
CHMAIR^,  ETWAC**,  NUTS^,  and  HICHEM^,  treat  the  beam  at  one  z-pcsition  in 
space.  It  is  clearly  impossible  to  put  this  much  detail  into  the  propagation 
code,  which  typically  requires  several  hundred  z-steps,  and  expect  to  lo  any 
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kind  of  parametric  study.  Therefore  some  of  the  details  of  a  complete  air 

chemistry  description  must  be  compromised.  The  VIPEH  model  attempts  to 

describe  the  major  processes  affecting  propagation  of  a  beam  similar  to  the 
ATA  beam  entering  the  atmosphere  in  a  way  which  Is  fast  enough  numerically  to 

be  useful. 

Th«  VIPER  chemistry  model  is  similar  in  structure  and  complexity  to  those 

used  in  the  SMPULSS'^  and  PHLAP^  propagation  codes.  This  model  has  been 

incorporated  into  other  NHL  propagation  codes  such  as  SIMMO,  SIM.M1 ,  and 

SARLAC  '  .  In  these  codes,  doth  chemistry  and  field  equations  involve  only 

derivatives  in  5.  The  VIPER  model  contains  monopole  and  dipole  rate  equations 

for  n^(r,;),  an  empirical  relation  for  calculating  T^  as  a  function  of  £/?, 

model  equations  for  rate  coefficients  and  momentum  transfer  frequency  v  as 

m 

functions  of  Tg,  and  the  usual  definition  a  -  n^e^/mv^.  The  original  VIPER 
model,  which  assumes  constant  v  ,  is  still  used  on  some  occasions. 

tU 

SEAM  PRODUCTION 

The  beams  of  interest  have  energies  in  the  MeV  range,  rthen  the  beam 
enters  a  gas  it  ionizes  the  molecules,  atoms  and  ions  present.  The  energetic 
secondaries  created  in  these  collisions,  in  turn,  produce  further  ionization 
and  energy  deposition.  If  the  mean  free  path  is  small,  the  secondaries 
cascade  down  in  energy  near  the  point  in  space  where  they  were  born.  This  is 
an  assumption  made  in  the  VIPER  model;  beam  ionization  and  the  associated 
cascade  is  iOcai  and  instantaneous.  This  is  a  good  approximation  in  760  torn 
air  and  is  probably  acceptable  down  to  -  100  torr’^.  3elow  that  nonlocal 
effects  begin  to  be  seen*^’^^.  Thus  beam  ionization  or  direct  production  is 
modeled  by 

dr. 

e 
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where  s  is  a  constant  direct  production  coefficient  (during  a  given  run),  P  is 
ratio  of  the  background  density  to  that  of  a  standard  atmosphere,  and  Jjj  is 
the  current  density  of  the  beam,  s  can  be  written  as 


s  •  (dc/dx)/(W)  (2) 

where  de/dx  is  the  stopping  power  in  a  standard  atmosphere  and  W  is  the 
average  energy  needed  to  create  an  electron-ion  pair.  Both  of  these 
quantities  depend  weakly  on  the  beam  energy  and  the  atmospheric  specie 
concentration.  In  molecular  and  atomic  nitrogen  and  oxygen  dc/dx  rises  from 
roughly  2200  eV/cm  for  1  MeV  beam  electrons  to  -  3000  eV/cm  for  100  MeV  beam 
electrons’^.  Though  this  variation  is  small  it  is  important  to  keep  in  mind 
when  comparing  predictions  for  different  accelerators.  Below  1  MeV  the  energy 
loss  increases  rapidly.  The  energy  per  ion  pair  is  around  35  eV  for  the  gases 
found  in  distrubed  air’^.  The  VIPER  value  used  for  ATA  (beam  energy  50  MeV) 
is: 

dn  /dt  Ll/cm'^  -  sec]  -  5.^  «  10  PJ^,^  [ amps/ cm  ]  (3) 

which  corresponds  to  (de/dx)/(DW)  -  36.4  [*/cm]. 

Direct  production  is  usually  the  major  conductivity  production  effect  and 
it  is  important  throughout  the  pulse. 

AVALANCHE  IONIZATION 

A  high  current,  fast  rising  beam  produces  large  electric  fields.  These 
fields  can  cause  the  plasma  electrons  to  avalanche,  producing  an  exponential 
growth  in  the  electron  density  and  conductivity  which  tries  to  short  out  the 
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field.  Large  radial  electric  fields  occur  in  the  r.cse  before  the  beam  is 
charged  neutralized.  A  large  inductive  axial  field  occurs  on  axis  if  the 
current  rise  time  is  short.  As  the  pressure  is  lowered  the  threshold  for 
avalanche  is  easier  to  achieve. 

It  has  been  shown  that  the  hollowing  instability  can  be  excited  by 

ft 

avalanche  .  The  VIPER  code,  because  it  uses  an  envelope  model,  cannot  study 
the  hollowing  instability. 

The  contribution  to  the  electron  density  due  to  avalanche  is  given  by 


dr.  /dt  -  a(E/P)n 
®  6 


(4) 


where  a  is  the  a-^alanche  coefficient  a.nd  E/P  is  the  magnitude  of  the  electric 
field  divided  by  the  pressure.  The  expression  for  a  is  fou.nd  in  Ali’^^ 


a/p 


[5x10*”^  ♦  0.21x10^(|)]  ex?(-273.8/(E/?))  for  50  <  E,  P  <  1 


20  V/ cm/ torn 


>.97x10^(|)’^2  exp(-359/(E/P)] 


for  E/P  >  120  v/om/torr 


(5) 

This  expression  agrees  with  the  Felsenthal  and  Proud  model  used  in  .RINGSEARER 

and  in  the  original  VIPER  chemistry  model  for  moderate  values  of  E/P  and  is 

xower  beyond  that.  The  reason  for  the  differences  are  given  in  (15). 

In  most  beams  w.here  avalanche  ionization  is  important,  hose  growth  tends 

to  be  more  violent  because  high  E/P  values  lead  to  a  reduction  in  v  the 

m  ’ 

momentum  transfer  frequency.  If  avalanche  is  strong  enough  to  excite  the 
hollowi.ng  instability^ ,  the  resulting  peaking  of  return  current  flow  or.  axis 
is  strongly  destabilizing  for  hose  as  well.  Thus,  most  VIPER  simulations  are 
performed  usi.ng  pa.-ameters  which  result  in  negligible  avalanche  ionization. 
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RECOMBINATION 

Dissociative  recombination  is  the  main  sink  for  electrons  in  disturbed 
air.  The  contribution  to  the  electron  density  is 

dn  /dt  -  “3  (6) 

e  e 


where  3  is  the  effective  recombination  coefficient. 

Recombination  is  not  important  in  the  nose  of  the  beam  where  the  electron 
density  is  low.  In  the  tail,  however,  it  balances  off  the  direct  production 
of  the  beam  and  the  electron  density  becomes  nearly  constant. 

Recombination  helps  to  stabilize  the  hose^"^,  because,  although  it  limits 
the  amount  of  conductivity,  it  broadens  the  radial  profile  of  the  conductivity 
and  the  plasma  current.  Plasma  current  flowing  outside  the  beam  has  a  strong 
stabilizing  effect  on  hose. 

Recombination  rates  are  strongly  dependent  on  the  type  of  molecular  ion 
present.  It  can  vary  by  twc  orders  of  magnitude  for  o.ne  volt  electrons.  In 
dry  air  the  predominant  ions  early  in  the  pulse  are  N^  and  These  react 
and  produce  many  cluster  ions  later  in  the  pulse.  The  cluster  ions  have 
faster  recombination  rates.  For  very  long  times  NO"^  dominates.  If  water 
vapor  is  present  the  effective  recombination  rate  is  large  because  of  the 
hydrated  clusters  which  form**.  Complicating  this  is  the  fact  that  many  of  the 
rates  are  not  well  verified  experimentally. 

The  recombination  rate  coefficient  used  in  VIPER  is 

3  -  0.66(N*)  -  Q.n&iol) 


(7) 


i 


wnere  3(Np  -  —  is  the  coefficient  for  M2*  and  3 (Op  -  —  is 


♦  13  '  ® 

that  for  Og  .  This  combination  was  chosen  to  fit  mere  sophisticated  air 


chemistry  modelling  results  for  fast  rising  beams  like  ATA  in  dry  air'^. 

In  summary,  recombination  must  be  included  if  the  bean  is  over  a  few 
nanoseconds  long.  It  is  very  dependent  on  the  beam  parameters  and  the 
compesitien  of  the  background.  Many  of  the  rates  are  net  well  measured.  A 
constant  rate  is  probably  acceptable,  but  if  temperature  dependence  is 
included,  the  rate  should  decrease  slowly  as  the  temperature  rises. 


ATTACKMSNT 

Oxygen  readily  attaches  low  energy  electrons.  Attachment  primarily 
occurs  in  the  wings  and  tail  of  a  beam.  It  is  a  three  body  process  so  it 
becomes  lass  important  s'!  the  density  decreases.  Attachment  is  not  included 
in  the  VIPER  model. 


lOLLISION  FRECUENCIES 


The  electron-neutral  collision  frequency,  comes  from  a  polynomial 


’it  to  rates  computed  for  a  Maxwellian  electron  distribution^^  using  a  recent 


compilation  of  the  momentum  transfer  cross  sections  for  nitrogen  and 
9  -ua  <•( ' 


5)f7gen' 


I  he  fit  is  given  by 


(  5.3966x10*^  *  1  .6629x10*'T^  -  2.5852x10*°  T*^  *  2.550Ux10*^T^  T  <  5  eV 
^  j  e  e  e  e 


.-3  ,2 


)  1 .1 153x10*'^  -  1. 7266x1 0*®' 


a. 3897x10 


■10  ,2 


>  5  eV 


(3) 


where  D  is  the  neutral  density  in  cm  Figure  ’  shows  the  polynomial  fit  for 
V  /D  and  the  computed  rates.  Since  T  <  2  eV  in  most  cases  of  interest,  the 
polynomial  fit  is  quits  accurate. 


The  electror.-ior.  collision  frequency,  is  given 
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ei 
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(9) 


The  electron-ion  contribution  is  negligible  in  the  nose.  For  long  beams 
the  ion  density  can  become  large  enough  so  that  the  electron-ion  term  can 
become  comparable  to,  or  even  surpass  the  electron-neutral  term.  In  the 
current  model  InA  is  held  fixed,  although  for  long  pulse  beams  its  variation 
should  be  considered. 

The  total  collision  frequency  is  simply  taken  to  be  the  sum  of  the 
electron-neutral  and  electron-ion  terms.  This  is  a  good  approximation  when 
one  of  them  dominates.  The  proper  weighting  requires  a  knowledge  of  the 
electron  distribution  function*'. 

The  conductivity  is  then  given  by 


e_ 
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en 


''ei^* 


(10) 


It  has  been  shown  that  inclusion  of  a  temperature  (or  E/P)  dependence  in 
the  collision  frequency  results  in  a  new  destabilizing  term  to  the  dipole 
conductivity.  When  compared  with  a  fixed  conductivity  model  an  ”E/P"  model 
usually  gives  more  growth  in  the  beam  head,  but  reduced  local  growth  rates  in 
the  beam  tail- ' . 


ELECTRON  TEMPERATURE 

It  is  assumed  that  the  plasma  sea  of  electrons  are  Maxwellian,  or  at 
least  can  be  characterized  by  an  electron  temperature.  There  are  many 


processes  that  ocntribute  to  the  determination  of  this  temperat'ure.  The  major 


The  electrons 


heating  processes  are  direct  team  heating  and  ohmic  deposition, 
which  enter  the  plasma  after  they  cascade  down  in  energy  from  their  creation 
in  the  beam  ionization  process  mentioned  above,  are  assumed  to  be  "bcrn"  with 
an  energy  of  about  7.5  eV .  This  source  of  energy  continues  as  long  as  the 
beam  is  on  and  is  the  dominant  heating  source  in  the  body.  Ohmic  heating, 
which  is  proportional  to  the  square  of  the  total  electric  field,  makes  a 
significant  contribution  in  t.he  nose.  The  electrons  are  cooled  by  their 
interactions  with  the  background  gas.  The  major  source  of  cooling  is  the 
excitation  of  the  vibrational  levels  of  nitrogen.  These  have  a  large  cross 
section  for  collisions  with  low  energy  electrons. 

Rather  than  try  to  solve  a  rate  equation  for  the  electron  temperature, 

the  <ll?ZR  model  uses  an  analytic  relation  between  the  local  electric  field  and 
the  temperature 


T 
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0.1 


,E.0.3 


*  T 

0 


E/?  in  vclts/cm/torr . 


(11  ) 


T..e  .  ^rst  term  comes  from  an  amalysis  of  ionisation  rates  bv  Ali^^.  ''t 
reflects  a  balance  between  ohmic  heating  and  excitation  cooling.  T^  is  a 
constant  which  models  the  beam  heating  component.  It  is  set  at  0.15  sV  for 
AiA  type  beams  at  fall  density.  For  lower  pressure  or  for  more  intense  beams 
a  higher  value,  typically  .30  -  .35  eV,  is  used. 

in  some  versions  a  relation  is  used  wnioh  was  found  from  a  fit  to  results 
from  the  air  chemistry  code  CHMAIH  11^- 


Tg  -  0.20  ‘  0.313  *  0.017a  I 


^  in  kv/cm/atm. 


(12) 
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The  Tg(E/P}  relationship  is  most  accurate  in  regions  of  the  beam  where 
ohmic  heating  is  the  dominant  heating  mechanism.  Thus,  it  is  well-suited  for 
treating  the  beam  head  and  body  where  most  hose  growth  usually  occurs.  It  is 
less  accurate  for  the  tail  region  of  long  pulses  where  collisional  beam 
deposition  can  be  the  dominant  heating  mechanism. 


SUMMARY  OF  MONOPOLE  AND  DIPOLE  EQUATIONS 

The  results  from  the  previous  sections  are  combined  to  give  rate 

* 

equations  for  the  monopcle  and  dipole  electron  densities  ng(r,c)  and  n  (r,;), 

Replacing  the  time  derivatives  in  (3),  (iJ),  and  (6)  with  c  gives 

oC 


c  -  sPJ.  *  a(E/P)n^  -  6(T^)n^ 
d^  b  e  e  e 


the  rates  s,  a,  and  S  are  given  in  (3),  (5),  and  (7)  respectively.  The 
conductivity  oCn^.T^)  is  given  in  (10),  using  the  collision  frequencies  in  (3) 
and  (9).  The  temperature  Tg  is  taken  from  the  E/P  relations  in  (11)  or  (12). 
The  monopole  conductivity  equations  are  used  to  give  oir,z)  in  VIPER,  SIMM0, 
and  SIMM1  and  a(r,;,e)  in  the  SARLAC  model. 

Linearized  hose  codes  such  as  VIPER  and  SIMMI  require  equations  for  the 
dipole  or  perturbed  conductivity.  Additional  terms  appea-  for  the  temperature 
dependent  rate  coefficients.  Using  to  denote  dipole  quantities,  we  have 


‘e  «  -  -  -  - 

—  -  sPJ  an  -  an  -  3r.  -  2Sn  n  , 

;  e  e  e  e  e’ 
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Since  the  coefficients  a,  3,  and  y  are  functions  of  T-(E/P),  the  dipole 

n  c 

A  A  A 

coefficients  a,  3,  and  are  given  by  appropriate  derivatives  of  these 
coefficients.  For  example, 


:  .  i«_  9  .  ii.  !Is  ; 

^  3T  e  3T  3E 
e  e 


with  the  dipole  electric  field  E  given  by 


E  E  +  E  E 
2  2  r  r 

(e2  >  ' 


(16) 


(17) 


The  actual  equations  used  in  VIPER  employ  a  set  of  dimensionless  units 
widely  used  in  the  beam  propagation  community.  Using  primes  to  denote 
dimensionless  quantities  and  scaling  to  the  nominal  beam  radius  a^,  we  have 
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The  dimensional  quantities  are  in  cgs  units,  and  -  e  /me  is  the  classical 
electron  radius.  The  dimensionless  forms  of  (10)  and  (iS'lS)  have  c  and  e^/m 
replaced  by  1.  The  detailed  diagnostics  output  from  VIPER  are  in  terms  of  the 
dimensionless  units  described  above. 


NUMERICAL  CONSIDERATIONS 

The  differential  equations  are  solved  by  a  simple  differencing.  If  the 
avalanching  becomes  important  an  integrating  factor  is  automatically  used. 

The  conductivity  and  field  solver  are  advanced  by  a  leapfrog  scheme.  If  the 
fields  change  too  rapidly  or  the  hose  motion  becomes  too  violent  an  automatic 
subgridding  of  the  timestep  is  performed.  Although  the  coefficients  are  giver, 
by  analytic  formulas,  we  usually  began  a  run  by  setting  up  tables  indexed  on 
E/P  a.nd  use  a  fast  table  lookup  routine  thereafter. 

BENCHMARK 

Table  1  shows  a  comparison  for  some  of  the  air  chemistry  quantities  for 
similar  runs  between  the  VIPER  model  discussed  above  and  the  radially  resolved 
CHMAIR  II  model^^.  CHMAIR  II  is  a  detailed  air  chemistry  code  containing 
many  species  and  reactions.  The  data  are  given  at  the  injection  point  for  a 
beam  similar  to  the  ATA  beam,  10  kA  and  50  MeV.  Figure  2  shows  the  electron 
temperature  plotted  versus  the  electric  field  at  several  radii  in  the  CHMAIR 
run.  The  VIPER  E/P  relationship  is  also  shown. 


Table  1.  Conparison  of  CHMAIH  II  and 


VIPER 

cliemlstry. 

Values  are 

given  on- 

'axis 

CItUIR 

VIPER 

CWAIR 

VI  PER 

CHKAIR 

VIPER 

Tia*(u) 

1 

1 

10 

10 

21  .6 

20 

Ej(v/sm) 

5350 

5827 

1J60 

1312 

1100 

376 

T,(*V) 

.58 

.66 

•  323 

.319 

.327 

.284 

9r, 

2U 

28.5 

14.3 

'3.3 

15 

15 

.«(<!' lasa) 
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3.0x10"’® 

2.9x10"’'^ 
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1.66x10"’® 

0 
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7.8x10"'^ 

7.7x10*'^ 

1.26x10"® 

1 .06x10"® 

9.3x10"' 

1 .1  xIO 

9( J'leaa) 

2.75 

2.36 

200 

233 

315 

234 
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CONCLUSIONS 


We  have  presented  an  air  chemistry  model  which  is  numerically  fast  and 
stable  enough  to  be  used  in  a  linearized  propagation  code.  The  model  is  tuned 
to  electron  beams  comparable  to  the  ATA  beam  entering  dry  air.  Processes 
determining  the  electron  density  are  direct  ionization,  avalanche  ionization, 
and  dissociative  recombination.  Temperature  dependent  electron-neutral  and 
electron-ion  collision  frequencies  are  included. 
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COLLISION  FR.  FOR  CONO 


Figure  1 . 


VIPER  polynomial  fit  compared  to  calculated  momentum  transfer 
frequency.  Tine  dasrted  curve  is  the  VIPER  fit,  the  dotted  curve 
is  a  weighted  sum  of  v  for  H2  and  O2,  and  the  solid  curve  is  an 
"exact"  calculation  of^v  using  the  integration  method  described 


VJIO) 
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vs.  E  for  CHMAI.R  run.  The  data  poi.nts  are  taken  at  r»0,  .25, 
and  1  cm  down  the  beam.  The  X's  are  what  VIPER  would  predict 


Figure  2. 
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